Critical properties of the three- and four-dimensional gauge glass 
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The gauge glass in dimensions three and four is studied using a variety of numerical methods 
in order to obtain accurate and reliable values for the critical parameters. Detailed comparisons 
are made of the sensitivity of the different techniques to corrections to finite-size scaling, which are 
generally the major source of systematic error in such measurements. For completeness we also 
present results in two dimensions. The variation of the critical exponents with space dimension is 
compared to results in Ising spin glasses. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
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INTRODUCTION 



Spin glasses have long been recognized to be the canon- 
ical examples of complex systems and so are of funda- 
mental interest for statistical physics and in particular 
the statistical physics of phase transitions. The critical 
behavior at standard continuous transitions in systems 
without disorder is now understood very precisely and 
in impressive detail thanks to the renormalization group 
(RG) approach. An important step forward towards the 
comprehension of the vast family of glassy transitions 
would be to arrive at a similar degree of understanding 
concerning the critical behavior at the spin-glass transi- 
tion. Unfortunately this task has proven to be extremely 
difficult. First of all, below the upper critical dimension 
(d = 6), the renormalization group theory that works so 
well for standard transitions has been found to present 
formidable technical difficulties^ and is of little guidance 
in lower dimensions. Second, on the purely practical side, 
numerical work is laborious and has tended to lead to im- 
precise values of ordering temperatures and critical expo- 
nents. Randomly disordered samples are microscopically 
inequivalent to each other and to obtain data that are 
truly representative of the global average, measurements 
must be made over large numbers of disorder realizations. 
Time scales are long and obtaining thermal equilibrium 
in large system sizes is hard. Finite-size scaling meth- 
ods are invariably used to estimate exponents but are 
carried out on a restricted range of sizes as the need for 
extensive computer time escalates when the system sizes 
increase. Because measurements are thus intrinsically re- 
stricted to small or moderate system sizes, it is essential 
to take into account possible artifacts due to corrections 
to finite-size scaling. Little is known about the magni- 
tude of these corrections a priori and nevertheless the 
data must be analyzed in such a way that existing cor- 
rections are identified and allowed for. 

Among spin glasses the Ising (ISG) systems have been 
by far the most studied numerically for obvious practi- 
cal reasons. Here we report results using a variety of 
numerical techniques on the gauge glass (GG) in space 
dimensions d = 3 and 4. This system is of interest in its 



own right because vector spin systems are generally much 
closer to experimental realizations of spin-glass ordering 
than are ISGs. Most laboratory spin glasses are made 
up of Heisenberg spins with 0(3)-symmetry and it has 
been suggested that chiral ordering plays an important 
role in spin freezingi^ Vortex glasses in type II super- 
conductors can be modeled by XY vector spins [0(2) 
symmetry], and GGs have already been studied exten- 
sively in this context ASiSii GGs are vector spin glasses 
which, due to symmetry arguments, do not show chiral 
ordering^ 

This work has two main aims. The first is to take this 
particular family of spin glasses as a case study in order to 
demonstrate that there exists a whole toolbox of numer- 
ical techniques available to identify spin-glass ordering 
transitions and to attempt to estimate as reliably and 
precisely as possible the related critical exponents and 
the associated corrections to scaling. Each of these tech- 
niques can have its advantages and disadvantages, and 
we evaluate the reliability of the different methods. Sec- 
ond, having values for the GG exponents in hand which 
are as precise as possible, it is of interest to follow the 
evolution of their values as a function of space dimen- 
sion in this particular family and to compare with values 
obtained in other families of spin glasses. 

We have checked for and analyzed corrections to finite- 
size scaling in the observables measured; the relative in- 
fluence of these corrections varies considerably from one 
observable to another. This can explain inconsistencies 
between the estimates of critical temperatures and crit- 
ical exponents in various spin glasses which have been 
reported. We conclude that the "current" defined later 
(related to the domain-wall stiffness at low temperatures) 
is little affected by corrections to finite-size scaling and 
thus provides reliable estimates of the ordering temper- 
ature T c . Measurements using the spin-glass suscep- 
tibility %(L, T), alone or together with nonequilibrium 
measurements, are also accurate and reliable. On the 
other hand in the systems with finite-temperature tran- 
sitions we have studied the intersection of the correlation 
length divided by system size, £l(T)/L, and find that 
the method is very sensitive to corrections to finite-size 
scaling even for sizes where these corrections have be- 
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come negligible for other observables. This method on 
its own would give imprecise or misleading estimates for 
the ordering temperature of the GG. Nevertheless it can 
relatively well pinpoint if T c > 0, or not, in general. 

In Sec. m we introduce the model studied. In Sec. IIIII 
we discuss the equilibration test used in the (parallel tem- 
pering) Monte Carlo method and describe the equilib- 
rium observables measured. Off-equilibrium Monte Carlo 
methods and observables are discussed in Sec. IIVI fol- 
lowed by results in d = 4 in Sec. M and d = 3 in Sec. I VII 
Some results for d — 2 are presented in Scc. lVIII We con- 
clude with a summary and a comparison of the techniques 
used to determine the critical parameters in Sec. lVIIll and 
with a discussion about the dimensional dependence of 
the critical exponents in Sec. IIXI Concluding remarks 
are presented in Sec. [X] 

II. MODEL 

The gauge glass is a canonical vector spin glass (see, 
for instance, Ref. 0) where XY spins on a [hyperjcubic 
lattice of size L interact through the Hamiltonian 

n = -jJ2 cos (^ " fa ~ A a) > (!) 

the sum ranging over nearest neighbors. The angles 4n 
represent the orientations of the XY spins and the Ay- 
are quenched random variables uniformly distributed be- 
tween [0, 27r] with the constraint that Aij — —Aji. J is 
conventionally set equal to 1. Periodic boundary condi- 
tions are applied. 

III. EQUILIBRIUM OBSERVABLES 

Equilibrium measurements are carried out with sam- 
ples fully thermalized using the exchange Monte Carlo 
(parallel tempering) technique&S We ensure equilibra- 
tion by checking that different observables do not change 
with the amount of Monte Carlo steps and measure by 
doubling the number of Monte Carlo steps. Once the 
last three measurements agree within error bars we are 
satisfied with the equilibration. 

For XY spin systems there is a choice to be made 
in the allowed single-spin acceptance angle for individ- 
ual updating steps. To optimize the updating procedure, 
the limiting angle is often chosen to be less than 2ir for 
an XY spinifl and linearly dependent on temperature. 
The numerical prefactor for the temperature-dependent 
window is chosen so that the acceptance ratios for the 
local Monte Carlo updates is ~ 0.4. As far as the final 
equilibrium parameters are concerned, this choice plays 
no role. However, for the noncquilibrium simulations to 
be introduced later on it is important to use the full 2ir 
acceptance angle. 

In Tabled we show N samp (number of samples), N sweep 
(total number of sweeps performed by each set of spins) , 



TABLE I: Parameters of the equilibrium simulations in four 
dimensions. iVsamp is the number of samples, A sweop is the to- 
tal number of Monte Carlo sweeps for each of the 2Nt copies 
(two replicas per temperature) for a single sample, and Nt 
is the number of temperatures used in the parallel tempering 
method. jVsam P (£.L) is the total number of of disorder real- 
izations used in the calculation of the two-point correlation 
length £l (defined below). The lowest temperature used is 
0.70, the highest 1.345. 



L 


N 

1 * samp 




TV 

1 ' sweep 


N T 


3 


5000 


1660 


2.0 x 10 4 


17 


4 


5000 


1250 


8.0 x 10 4 


17 


5 


5000 


1000 


4.0 x 10 5 


17 



TABLE II: Parameters of the equilibrium simulations in three 
dimensions. The lowest temperature simulated is 0.05, the 
highest 0.947. The different quantities are explained in the 
caption of Table Q 



L 


iVsamp 


Asamp(Ci) 


^Vsweep 


N T 


3 


10000 


2660 


6.0 x 10 3 


53 


4 


10000 


2000 


2.0 x 10 4 


53 


5 


10000 


1600 


6.0 x 10 4 


53 


6 


5000 


1330 


2.0 x 10 5 


53 


8 


2000 


1000 


1.2 x 10 6 


53 



and Nt (number of temperature values) , used in the sim- 
ulations in four dimensions. Table ITT1 has the correspond- 
ing values for the simulations in three dimensions. The 
parameters for the simulations in two dimensions are pre- 
sented in Table IIIII (see also Ref. . 

A primary observable for a spin glass is the equilibrium 
spin-glass susceptibility at finite system size L, x(L,T). 
The susceptibility is defined asm 

X = Niiq 2 )}^ , (2) 

where q is the spin-glass order parameter: 

1 N 

i=l 



TABLE III: Parameters of the equilibrium simulations in two 
dimensions. The lowest temperature used is 0.13, the highest 
1.058. For L — 24 the lowest temperature studied is 0.20. 
The different quantities are explained in the caption of Fable 



HI 


L 


-A/samp 


A r samp(£i) 


^Vsweep 


Nt 


4 


10400 


4000 


8.0 x 10 4 


30 


6 


10150 


2660 


8.0 x 10 4 


30 


8 


8495 


2000 


2.0 x 10 5 


30 


12 


6890 


1330 


8.0 x 10 5 


30 


16 


2500 


1000 


2.0 x 10 6 


30 


24 


2166 




2.0 x 10 6 


24 
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Here, a and (5 are two replicas of the system with the 
same disorder and [• • -] av represents a disorder average, 
whereas (• ■ •) represents a thermal average; N is the num- 
ber of spins in the system. The standard finite-size scal- 
ing form for the equilibrium spin-glass susceptibility is 

x(L) = L*-"C[L 1 ' v (T-T c )]. (4) 

The function C tends to a constant at the critical tem- 
perature so that at T c 

X (L) cx L 2 -" , (5) 

where rj and v are the usual critical exponents. Follow- 
ing the RG approach, the spin-glass susceptibility, like 
the other observables, will be modified at small L by a 
correction to scaling factor ~ [1 + AL~ U + •••], where uj 
is the correction exponent arising from the leading irrele- 
vant operator in the RG and A is a constant. However, in 
general there are also "lattice artifact" correction termsi^ 
giving at T c 

X (L) cx L 2 -" + B + • • • , (6) 

where B has no relation to the RG correction. This leads 
to an effective leading correction factor [1+B'L~( 2 ~' ? )] . If 
u> > (2 — 77) corrections will be dominated by the leading 
"analytic" term with an effective correction exponent (2— 

^ 12.13.14 

Spin-glass susceptibility measurements can be used di- 
rectly to estimate T c . As x{L) in absence of corrections 
to scaling increases as L 2 ~ v at the ordering temperature, 
a log-log plolpi^ of x{L) against L is linear at T c . At 
higher temperatures the plot will curve downward and at 
lower temperatures it will curve upward. The crossover 
from negative to positive curvature at large L should 
give a precise measure of T c . Alternatively one can plot 
\n[x(nL) / x(L)} as functions of T for different L, and, for 
instance, n = 2. The curves will intersect at T c . Direct 
scaling of x(L, T) according to Eq. |0J on the contrary 
provides a very poor indication for T c £ In the presence of 
corrections there will be an additional curvature at small 
L which should be essentially temperature independent 
near T c . If data on a reasonably wide range of L are avail- 
able it is possible to identify the leading correction term 
and to estimate T c using data at L values large enough 
for the correction term to be negligible. 

The Binder ratic 16 for the GG is defined by 17 



Finite-size scaling predicts 

g(T,L) = G[L 1 /"(T-T c )] . (8) 

This dimensionless observable should be independent of 
L at T c , except for corrections to scaling. The Binder 
ratio is a bona fide "work horse" widely used to obtain an 
estimate of the critical temperature of statistical systems. 



It has frequently been used in ISG studies to estimate 
ordering temperatures. However, even in ISGs the Binder 
crossing point method is not efficient for determining T c 
precisely, at least in dimension three, because the g(T, L) 
curves tend to lie very close to each other so the estimate 
of the crossing point is very sensitive to corrections to 
scaling. 18 In the particular case of the GG, the method 
is inoperable because the Binder curves do not intersect, 
at least for the range of sizes that have been studied here. 

For vector systems it is possible to measure a domain- 
wall-stiffness-like parameter, the "current" which 
is the rate of change of the free energy with respect to a 
twist angle at the boundaries. For the gauge glass: 

I N 

I{L) = — ^ sin(0j - cj) i+Sl - A l . (9) 
i=i 

In this case, the twist is applied along the x direction. 
As [(J(L))] av = 0, we actually calculate the root-mean- 
square current I rms between two replicas: 

7 rms = y/[(I(L) a )(I(L)p)U . (10) 

The root-mean-square currents for different L cross at T c 
and splay out at lower temperatures. They have a finite- 
size scaling form for a finite-temperature transition 7 

/ = /[L 1 /-(r-T c )] . (11) 

The intersection of the J rms (L) curves gives a clear indi- 
cation of the value of T c . Corrections to finite-size scal- 
ing must be allowed for, but in practice these are weak. 
The main drawback of this method is the fact that in 
current measurements intrinsic sample-to-sample fluctu- 
ations are very strong so that data must be taken on 
a large number of independent disorder realizations in 
thermal equilibrium. 

A further important observable is the two-point corre- 
lation length. In an infinite sample the correlation func- 
tion 

G(r ij ) = [{S i -S j ) 2 ] a;v (12) 

is of the form 

G( rij )o,r^ d - 2+ ^e- r ^ T \ (13) 

where i and j represent the position of the magnetic mo- 
ments. £,(T) is the correlation length, which diverges 
as (T — T c )~ v . In some systems G(r) has been directly 
recorded; when T c is known, the measurement of G(r) 
at T c on reasonably large samples provides a very di- 
rect measurement of 77. (Measurements have also been 
made as a function of anneal time for temperatures lower 
than T c in the Edwards Anderson ISG. 19 ) Under periodic 
boundary conditions there must clearly be a cutoff for 
G(r) at rij = L/2, and even before this cutoff the behav- 
ior is modified as G{rij) must become r-independent at 
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the cell boundary so as to be compatible with the peri- 
odic boundary conditions. A size-dependent correlation 
length £l can be defined through 20 ' 21 - 22 - 23 



Sl = 



1 



2sin(|k min |/2) 



C(k m in) 



1/2 



(14) 



where G(0) and G(k m ; n ) are Fourier transforms of 
the spatial correlation function G(r) 2 . 4 , and k m j n = 
(27T/.L, 0, 0) is the smallest nonzero wave vector. This 
second-moment correlation length is in fact an observ- 
able having the dimension of length, which becomes equal 
to the correlation length in the limit L £. It is re- 
ferred to as the "size-dependent correlation length" even 
in the limit of T = T c and below, where the infinite- 
size correlation length £ has diverged. The ratio £l /L at 
T c should be L-independen1*ii as the form of the whole 
function G(r) scales appropriately with L. Curves for the 
ratio of the finite-size correlation length to sample length 
£l (T) / L for different L cross at T c and then splay out at 
lower T, 11 ! 20 ! 21 ! 25 i.e., 



U/L = X[L 1 '"(T-T e )] 



(15) 



This makes this observable very attractive for estimat- 
ing T c . Unfortunately at small and moderate L, the ratio 
£l/L can be very susceptible to corrections to scaling. In 
addition, the definition of £l through Eq. (|14|) is a con- 
vention and other definitions having the same limiting 
form at infinite L are equally plausible^ At small and 
moderate L there are correction terms ~ L~ 2 from the 
subleading term in the sine factor in the conventional 
definition. There may be a "lattice artifact" correction 
tern: 12 of order 1/L, in addition to the "true" correction 
term in L~ w arising from the leading irrelevant operator. 
Physically, G(0) is the sum over the spins within a box 
~ L d correlated to a central spin and so is equal to the 
equilibrium x{L)- G(2tt/L) contains a positive term from 
spins close to the central spin together with a term from 
some of the spins further from the central spin which is 
negative because of the cosine factor. A slight change 
in the form of the function G(r) with L will have little 
effect on the L-dependence of G(0) while it can modify 
the Fourier transform G(2n/L) much more drastically. 
Hence one can expect that at low and moderate L, £l/L 
will be much more sensitive to corrections to finite-size 
scaling than x(-^) will. Even for the canonical ferromag- 
netic Ising model in two dimensions, the corrections in 
£l/ f L (or to the variant second moment correlation length 
£' L /L with a slightly different definition 1 ^) at criticality 
are strong for L < 20, with different correction terms 
contributing^ As a consequence, correlation- length data 
should as a general rule be treated with caution unless 
detailed results exist over a wide range of L from which 
the various correction contributions can be evaluated. 

Taken over the entire temperature range, ratios such 
as &l{T)/£ l {T) or x(2i, T)/ X (L,T) should be universal 
scaling functions of £l(T) / Lm&H there are corrections to 



scaling, scaling curves for different L will not be identi- 
cal; this provides an operational method for checking for 
corrections to scaling in £l/L not only at or near T c , but 
over the whole range of temperatures. 

Finally, we also record the equilibrium energy per spin. 
As is well known, the energy of glassy systems varies 
smoothly through the ordering temperature, but the size 
dependence of the energy has a behavior linked to the 
ordering temperature. This is discussed elsewhere. 26 



IV. OUT OF EQUILIBRIUM METHODS 

Complementary out-of-equilibrium measurements are 
also carried out. Large samples of size L are initialized 
in a disordered (infinite temperature) state; they are then 
held at a bath temperature T for an anneal time t w . The 
spin-glass susceptibility x(^^w) is monitored as a func- 
tion of t. At T c dynamic scaling rules hold; the spin-glass 
susceptibility increases proportional to i( 2-1 )/ 2 ), as clus- 
ters of correlated spins build up over timei2£ This equa- 
tion is not valid at very short time scales as \ is identi- 
cally equal to 1 for a sample of any size and dimension in 
the totally disordered state, so the scaling equation will 
only set in after a crossover from this initial zero-time 
value. Even after this crossover, there can be corrections 
to scaling in the nonequilibrium susceptibility data, when 
only small clusters have been built up. Here the correc- 
tions appear as corrections to finite-time scaling but are 
physically equivalent to corrections to finite-size scaling. 
Including the correction term, 



x(t) = A<( 2 -"/^(i + Br w ' z + ■■■) 



(16) 



where w is a dynamic correction to scaling exponent anal- 
ogous to uj In principle w and u> have no fundamental 
reason to be exactly equal but can be expected to be 
similar. Finally, if the measurement is continued long 
enough, the long-time susceptibility will cross over to a 
saturation value which is just the equilibrium suscepti- 
bility x(L,T) for the measurement temperature T and 
the sample size L used. 

It is clear from a comparison of the equations for the 
equilibrium susceptibility and the dynamic susceptibility 
at T c that for all intermediate times between the very 
short time limit and the L-dependent long-time limit, the 
measurement after a time t w on a large sample should 
be equivalent to the measurement at equilibrium on a 
sample of size 



L* = At 1/z 



(17) 



with A a constant and z the dynamic scaling exponent. 
x(L) and x(L*) can be displayed on the same graph; with 
a judicious choice of parameters A(T) and z{T) the set of 
x{L*) scale onto the set of equilibrium data x(£). This 
provides us with a direct method for estimating z(T). Al- 
ternative techniques which have been used for estimating 
the value of z(T) independently of the other exponents 
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are the monitoring of the time variation of the Binder pa- 
rameter for different sample sizes^i or the scaling of the 
time-dependent spin-glass susceptibility x{L,t) for dif- 
ferent sizes The present method is more convenient 
than the Binder parameter measurements as the latter 
are intrinsically noisy. 

We have carried out dynamical measurements over a 
range of temperatures and not just at the putative order- 
ing temperature. If we assume that there is an effective 
dynamic scaling exponent z(T) at each temperature T 
and not only at T c , then at each T we should be able to 
translate x(t) data into %(£*) data by a suitable choice 
of A(T) and z(T). If dynamic scaling continues to hold 
at other T, then once again it should be possible to make 
the x{L*) data scale onto the x(£) data. The assumption 
of an effective z(T) has been made before in a different 
context, but applied only at temperatures below T c (see, 
for instance. Refs. and 

Finally, after a long anneal time t w the ultimate config- 
uration {Si(£ w )} is registered. The updating procedure 
is then pursued for a further time t considerably less than 
t w and the autocorrelation function decay 

1 - 

g(t,t w ) = ^2Z(Si(t w +t)-S i (t w )> (18) 

is monitored starting from what can now be considered 
a quasi-equilibrium configuration. In Eq. I|18|) repre- 
sents the vector spins in the plane. In the context of 
Ising spin glasses it has been showni2i2ii22 that the ini- 
tial decay behavior of q(t, t w ) is not dependent on the 
sample having achieved perfect thermal equilibrium; as 
long as the time scale over which q(t, t w ) is monitored is 
much shorter than the annealing time i w , then the form 
of q(t,t w ) is characteristic of the infinite-size limit initial 
relaxation in equilibrium. Again at T c , for a large sam- 
ple in equilibrium, the autocorrelation function q(t,t w ) 
decays algebraically as q(t,t w ) ~ t~ x with^ 3 . 

x = {d-2 + ri)/2z {T = T C ). (19) 

Measurements in Ising spin glassesi^i 3 ^* 3 ^ 3 ^ have shown 
that at all temperatures below T c the q(t, f w ) data can 
be accurately fitted by this power-law form of decay with 
a temperature-dependent effective exponent x(T). For 
temperatures higher than T c , the decay of q(t, t w ) has the 
same initial power-law form which is now multiplied by 
a further factor Q[t/r(T)], where Q is a scaling function 
and t(T) is a relaxation time diverging as (T — T c )~ zv 
when T c is approached from above^ In ISG measure- 
ments q(t,t w ) shows very little short-time corrections to 
scaling; at T = T c and below, q(t 7 t w ) for large well- 
annealed samples follows a strict power-law decay beyond 
very few Monte Carlo steps after "zero" time. 

Nonequilibrium measurements can be used in combina- 
tion with the equilibrium susceptibility measurements to 
obtain an estimate of T c from consistency arguments. 35 
The effective dynamical exponent z(T) is measured from 



TABLE IV: Parameters of the off-equilibrium simulations as 
a function of space dimension d. L is the size of the system 
used, t the equilibration time, and t w is the "waiting time" to 
calculate q(t, £ w ), Eq. 1181 . JV samp is the number of samples 
used for the disorder average. 



d 


L 


■^Vj samp 


t 




2 


64 


500 


1.638 x 10 4 


4.069 x 10 3 


3 


16 


200 


8.192 x 10 3 


8.192 x 10 3 


4 


10 


200 


1.310 x 10 s 


4.096 x 10 3 



the comparison of equilibrium and nonequilibrium sus- 
ceptibility measurements. The effective exponent rj(T) 
is measured from the size-dependent equilibrium sus- 
ceptibility measurements through x(L,T)/L 2 — L~ V ( T \ 
Then, at T c , the relaxation decay exponent x must be 
equal to (d — 2 + rj)/2z. The consistency between the 
directly measured x(T) and the other exponents gives an 
efficient criterion for determining T c . This technique has 
a number of advantages. x(T) and z{T) can be mea- 
sured on almost arbitrarily large samples and so can be 
rendered free of corrections to scaling. In addition, they 
can be measured without the need to achieve complete 
thermal equilibrium. The most laborious measurement 
is the size-dependence of xCO m equilibrium; however, 
x{L, T) measurements are generally only subject to weak 
corrections to finite-size scaling so a limited range of L 
can be adequate. 

The parameters of the off-equilibrium simulations us- 
ing simple Monte Carlo updates are listed in Tabic Hvl 



V. FOUR DIMENSIONS 

In this section we compare different observables as 
well as critical exponents and transition temperatures 
derived from finite-size scaling arguments for the four- 
dimensional gauge glass. 



A. Root-mean-square current 

Current measurements are made for sizes L — 3 - 
5. The details of the simulation are shown in Table [I] 
The intersection of the J rms curves shown in Fig. Ogives 
an estimate T c = 0.890 ± 0.015. This is significantly 
lower than the value of 0.96 ± 0.01 reported by Reger 
and Young* 5 , although their data are almost compatible 
with the present intersection point within the statistical 
error bars. Corrections to finite-size scaling appear to 
be minimal for these measurements. However, current 
measurements are intrinsically noisy with strong samplc- 
to-sample variations at equilibrium, so it is imperative to 
average over large numbers of samples. In Fig.|21we show 
a scaling plot of the data in Fig. ^ according to Eq. Qlip. 
The data scale well for the (modest) range of sizes and 
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FIG. 1: (Color online) Data for the the root-mean-square 
current I IU1S in four dimensions. The data show a crossing at 
T c » 0.89. 
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FIG. 3: (Color online) Data for ln(x/£ 2 ) vs ln(L) for different 
temperatures in four dimensions. At T c one expects \ ~ 
[i-n -pjje data show a change in curvature while scanning 
through T c w 0.89. 
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FIG. 2: (Color online) Scaling of the root-mean-square cur- 
rent 7 rma in four dimensions according to Eq. 11 U . The 
data scale well for the range of sizes shown and we estimate 
1/u » 1.42 with T c w 0.89. 



we estimate 1/v = 1.42 ± 0.03 together with the above- 
mentioned estimate of T c . The previously quoted error 
is estimated by varying the scaling parameters until the 
data do not collapse well. This method is also used in all 
subsequent estimates of error bars of scaling exponents 
and critical temperatures derived from scaling plots. 



B. Equilibrium susceptibility 

The susceptibility data at equilibrium are obtained 
from the same data set in Sec. IV Al and shown in Fig. 
The data are plotted in the form of curves for ln[x(L) / L 2 ] 
against InL at different temperatures. As discussed 
above, if corrections to scaling are negligible, this form of 
plot gives a straight line of slope —rj at T c , together with 
curves which bend downward for T > T c and upward for 
T < T c . For the four-dimensional GG, corrections are 
very weak because even including values for L = 2 the 
log-log plot of the data at the temperature closest to T c 
as estimated above follows the straight line behavior very 
closely. 

In order to obtain an independent estimate for T c , we 
have plotted the x 2 -deviation from a straight line fit to 
the data over a range of temperatures around T = 0.89, 
together with the curvature for three-parameter fits, 
Figs. 01 and |SJ The results indicate straight-line behavior 
for T = 0.895 ± 0.015 which is an independent measure- 
ment of T c . In Fig. |S1 we also show the effective slope 
—rj(T) from the three-parameter fit. The slope of the 
straight line at T c gives an estimate of the critical expo- 
nent r] = —0.74 ± 0.03. As far as we are aware of, this 
is the first published estimate for this parameter for the 
GG in four dimensions. 

With values for T c and rj in hand, the whole 
X(L,T) data set can be plotted in a standard man- 
ner: X {L,T)/L 2 -v as a function of - T c ), ad- 
justing v to obtain optimal scaling. The results lead to 
1/V = 1.42 ± 0.03, as shown in Fig. and are in agree- 
ment with an estimate from the scaling of the currents 
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FIG. 4: (Color online) x 2- deviation from a straight-line fit 
for the susceptibility data in four dimensions shown in Fig. [3] 
The data show a minimum at T ~ 0.89, indicating that the 
optimal fit happens around T — T c (marked by an arrow). 



presented in Fig. [21 



C. Correlation length 



FIG. 5: (Color online) Data for the effective exponent r/ D s (T) 
as a function of temperature in four dimensions. The vertical 
dashed line marks our estimate of T c from current and suscep- 
tibility measurements, T c = 0.89. In addition, rf is shown, 
which is an effective exponent derived from off-equilibrium 
calculations (described in detail in Sec. llVII . One expects J? e ff 
and rf to cross at T c , which is the case in our data within 
error bars. Data for the curvature of the susceptibility shown 
in Fig. from a second-order polynomial fit are also dis- 
played. The data cross zero curvature (horizontal dashed line) 
at T w 0.90, a value slightly higher than the other estimates. 



The data for the ratio £l(T)/L for different sizes are 
shown in Fig. Although the data appear to cleanly 
intersect at one point, it can be seen that the intersection 
temperatures for these sizes are around T « 1.0, and are 
only approaching the true ordering temperature, T c — 
0.89, very slowly with increasing L. The comparison with 
the behavior observed for xO^O or the current is striking. 
For each of these two parameters the critical behavior at 
T c is virtually correction-free, while for the same range 
of small values of L, the £,l/L intersections give a very 
poor indication of the true ordering temperature. There 
are clearly strong corrections to scaling at small L which 
appear to be intrinsic to this form of measurement. 



D. Nonequilibrium susceptibility 

Large samples (L = 10 in four dimensions) are ini- 
tialized in random (infinite temperature) configurations. 
They are then put in contact with a heat bath at a tem- 
perature T, and the spin-glass susceptibility is recorded 
as a function of annealing time. At T c the spin-glass 
susceptibility will increase as f( 2-T ?)/ z until the suscepti- 
bility arrives at the equilibrium value for that size L. We 
present the data in terms of the scaled time-dependent 
effective length L* — At x / Z at each temperature. When 
A and z are suitably chosen, the equilibrium suscepti- 
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FIG. 6: (Color online) Scaling plot according to Eq. Q of 
the susceptibility data presented in Fig.|3]in four dimensions. 
The data scale well for T c w 0.89, 1/u « 1.42, and r\ ft* -0.74. 
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8 




10 3 



0.7 0. 



0.9 



1 
T 



1.1 1.2 1.3 



FIG. 7: (Color online) Correlation length £l divided by L for 
different system sizes in four dimensions. The data cross at 
T ~ 1.0 indicating that there are strong corrections to scaling. 
The inset shows a zoom of the data in the main panel focusing 
on the crossing point around T ~ 1. 



bility x(X) and the nonequilibrium susceptibility x(L*) 
before saturation scale well together. An example for a 
temperature close to T c , T = 0.894, is shown in Fig. |S] 
From the data, the dynamical scaling exponent z(T) can 
be determined accurately, see Fig. |5J The behavior of 
the effective dynamical exponent z(T) shows no appar- 
ent special behavior at the critical temperature. Effective 
temperature-dependent values z(T) have been reported 
in ISGs for T < T c (sec, for instance, Ref . l30|) . We can 
conclude from the present data for the GG that z(T) is a 
well-defined parameter for a whole range of temperatures 
including T > T c . At T = T c we estimate z = 4.50±0.05. 



E. Autocorrelation function decay 

In the GG we have found unexpectedly that q(t) only 
assumes a pure power-law behavior after a relatively long 
time, of the order of 100 Monte Carlo steps, hereafter re- 
ferred to as MCS, (as compared with ~ 5 MCS in ISGs), 
see Fig. ED For earlier times the decay is affected by 
short-time corrections. In practice this means that the 
measurements of x(T) as defined above are less precise 
than in ISGs as they are limited at short times by the 
correction term and at long times by the condition that 
the maximum t w should be much less than t. 

We compare, for a set of temperatures around T c , 
the value of rj(T) directly obtained from the equilib- 
rium susceptibility measurements with a value calcu- 
lated indirectly from the x(T) and z(T) data, rf(T) = 
(d-2)-2z(T)x(T). We expect rj(T) = rj*(T) at T = T c . 
This can be seen in Fig. [3] where the vertical dashed 
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FIG. 8: (Color online) Spin-glass susceptibility from equi- 
librium measurements x(L) at T = 0.894 in four dimen- 
sions. x{L") is the susceptibility determined from the off- 
equilibrium simulations with L* = At 1/Iz , A = 1.22, and 
z ~ 4.5. By suitably choosing L* the data for x(L) and 
X(L*) fall on a straight line. This allows us to determine the 
dynamical critical effective exponent z(T) as a function of 
temperature. When L*(t) approaches the sample size (here 
L = 10), x(L*) necessarily saturates. The dashed line is a 
guide to the eye. 
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FIG. 9: (Color online) Dynamical critical exponent z as a 
function of 1/T in four dimensions as determined from the 
procedure described in the text and in the caption of Fig. |H| 
The data are consistent with z(T = T c ) = 4.50 ± 0.05. Note 
that the error is estimated "by eye" : the data are varied until 
the two expressions for the susceptibility differ noticeably thus 
allowing us to give an upper bound for the error bars. The 
dashed line is a guide to the eye. 
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FIG. 10: (Color online) Data for the autocorrelation function 
q(t) in four dimensions as a function of time (measured in 
Monte Carlo steps). The dashed line has slope —x = —0.147 
from which we can determine rf , a "dynamical" effective ex- 
ponent r\. Details are described in the text. 



line marks our estimate of T c . The two curves cross at 
T = 0.90±0.01, r] = -0.74±0.02, which again defines the 
ordering temperature and the critical exponent r\. As for 
this particular system T c is already reliably established 
from the previously mentioned methods, the comparison 
provides a critical benchmark test for the nonequilibrium 
technique. The good agreement with the standard meth- 
ods, such as susceptibility scaling and the crossing of 
root-mean-square currents in this case shows that this 
method is reliable, implying that it can be used for other 
systems. 



VI. THREE DIMENSIONS 

A number of estimates have already been given of the 
ordering temperature and critical exponents of the GG in 
dimension threeA In this section we present estimates for 
the critical temperature and exponents as derived from 
our simulations. In comparison to the four-dimensional 
data presented in Sec. [V] we find strong corrections to 
scaling in three dimensions. In what follows the different 
results for different obscrvables are presented. 



A. Root-mean-square current 

Olson and Youn g& estimated that T c = 0.47 ± 0.03 
from the intersection of the root-mean-square current in- 
duced by an infinitesimal twist along the boundaries. The 
present current data are shown in Fig. 1111 Our results 
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FIG. 11: (Color online) Root-mean-square current 7 rms as a 
function of temperature for the three-dimensional GG. The 
data show a crossing at T ~ 0.46. Note that the data for 
L — 3 show strong corrections to scaling. 



agree with the data of Ref. within error bars. A con- 
sistent intersection point for all the data taken together 
occurs at T — 0.46 ± 0.01, with only the data for L = 3 
lying slightly below the region of intersection. We con- 
clude that this temperature indeed represents the cor- 
rect ordering temperature for the system, corrections to 
finite-size scaling being small for L > 4. The precision is 
limited mainly by the statistical error bars, as the curves 
for different sizes lie rather close together. A scaling plot 
of the root-mean-square current according to Eq. Ijll|l is 
shown in Fig.[T21for l/v = 0.72 and T c = 0.46. 



B. Equilibrium susceptibility 

A log-log plot of the equilibrium susceptibility divided 
by L 2 as function of L for different temperatures is shown 
in Fig. In absence of corrections to scaling this plot 
should be straight at T c . We see that for T > 0.46 the 
data curve down, indicative that we are above the critical 
temperature. For T < 0.46 the data show a slight upward 
curvature^ in agreement with T < T c . Using Eq. (J5j) at 
T = T c we obtain -q = —0.47 ± 0.02 and a correction 
factor - (1 - AL-") with lu ps 2.5 and A = 0.6. For 
temperatures below T c a higher effective rj{T) is obtained 
at each temperature and the correction factor appears 
to be temperature-independent. The correction term is 
almost negligible for L > 4; if we make straight line fits 
to ln[x(T)/L 2 } against InL for L = 4 to 8, we find x 2 
values which increase sharply above T ps 0.40, and which 
rise slowly for temperatures below T ps 0.40, see Fig. 1141 

In both Ising spin glasses and GGs the first term in 
the RG e-expansion for the leading irrelevant operator 
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FIG. 12: (Color online) Scaling of the root-mean-square cur- 
rent 7rms in three dimensions according to Eq. 1111 , We see 
acceptable scaling of the data around T = 0.46. Deviations at 
higher T are presumably due to corrections to scaling. This 
plot is for l/v = 0.72 and T c = 0.46. 
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FIG. 13: (Color online) Data for \n(x/L 2 ) vs In L for different 
temperatures in three dimensions. The data show a change 
in curvature in the log-log plot while scanning through the T c 
estimate from the root-mean-square currents, 0.46. At T c we 
expect x ~ L 2 ~ v . 



is u(d) ~ (6 — <i)2LiS implying that ui will be high in 
dimension three. High-temperature series expansion data 
show that u> is greater than 3 in dimension three for the 
ISG2&. If we make the plausible assumption that this is 
the case also for the GG, the leading correction term will 
be due to the lattice artifact so the correction factor is 
[1 + AL~( 2 -'')]. Fits to all 3D GG x(L,T) data from L = 
2 to 8 using this correction factor give a clear minimum 
in x 2 as a function of temperature at T = 0.45 ± 0.02 
(inset to Fig. I14|) . which we can identify with T c . The 
effective exponent rj c g(T) is shown in Fig. 1161 

An overall scaling plot of x{L, T)/ L 2 ^ v as a function of 
Li/v(T _ jy f or L > 4 is shown in Fig.UHl We estimate 
l/v = 0.72 ±0.02. 



C. Binder parameter 

In spin glasses the Binder parameter g defined in 
Eq. Q is independent of system size at the transition 
temperature as it is proportional to a function which only 
depends on h 1 ^ (T — T c ). Consequently different lines at 
different temperatures cross at T c . In vector systems this 
is not the case: data for different system sizes splay for 
T > T c but not for T < T c , as can be seen in Fig. IT71 
From the data one can, at best, estimate T c roughly be- 
cause the data do not cross. The same behavior is found 
in the three-dimensional XY spin glass^ This has been 
ascribed to the ordering being chiral in this model. As 
no chirality can be defined for the gauge glass, this is not 
the case here and the behavior of g for the GG remains 
to be understood. 
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FIG. 14: (Color online) x 2 -deviation from a straight-line fit 
for the susceptibility data in three dimensions. The main 
panel shows data for a fit with no corrections to scaling due to 
lattice artifacts, whereas in the inset we show the data where 
corrections to scaling are included in the fits of ln(x/L 2 ). The 
data in the inset show a minimum at T w 0.45, in agreement 
with data from Ref. and with the root-mean-square current 
estimate from Sec. I VI A1 
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FIG. 15: (Color online) Scaling plot according to Eq. IgJ of 
the susceptibility data presented in Fig. 1131 in three dimen- 
sions. Optimal scaling is obtained for T c ~ 0.46, 1/u « 0.72, 
and tj » -0.47. 




FIG. 17: (Color online) Data for the Binder ratio g as a func- 
tion of temperature for several system sizes. We see that the 
data do not cross for T — T c « 0.46 and also do not splay for 
T smaller than T c . The inset zooms into the region T < 0.6. 
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FIG. 16: (Color online) Data for the effective exponent »7eff(T) 
as a function of temperature in three dimensions. Also plot- 
ted is rf as estimated from off-equilibrium simulations. Both 
exponents do not cross cleanly making it difficult to determine 
T c . The vertical dashed line marks our estimate for T c , 0.46. 
Also shown is the curvature of the susceptibility data (see 
Fig. IIMI from a second-order polynomial fit. The data cross 
zero curvature (horizontal dashed line) at T = 0.44 ± 0.01, a 
value slightly lower than other estimates of T c . 



D. Correlation length 

The data for the ratio £l{T)/L are shown in Fig. ED 
As in four dimensions there are strong finite-size correc- 
tion effects; in addition, the curves for different L lie close 
together. The intersections between curves for different 
L change with temperature in such a way that from these 
data alone it would be very hard to identify T c to better 
than T c = 0.50 ± 0.05. 



E. Nonequilibrium susceptibility 

In just the same way as in four dimensions the 
nonequilibrium susceptibility is recorded for large sys- 
tem sizes, L — 16. The effective dynamical exponent 
z(T) is estimated from the comparison of the time- 
dependent nonequilibrium susceptibility and the size- 
dependent equilibrium susceptibility, Fig. El an d the ef- 
fective dynamical critical exponent z(T) is shown as a 
function of temperature in Fig. [201 At T = 0.46 we ob- 
tain z(T) = 4.7 ±0.1, in agreement with the results from 
Ref. but with smaller error bars. 



F. Autocorrelation function decay 

Corrections to finite time scaling in three dimensions 
extend to times of the order of 100 MCS, limiting the 
precision of the measurement of the exponent x(T), as in 
four dimensions. For this particular system the directly 
measured effective exponent rj(T) and the indirectly es- 
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FIG. 18: (Color online) Correlation length £l divided by L 
for different system sizes in three dimensions. The data cross 
at T ~ 0.5 but only splay slightly making it difficult to give 
a precise estimate of the transition temperature. 




FIG. 20: (Color online) Dynamical critical exponent 2 as a 
function of 1/T in three dimensions. For details see the text. 
The data are consistent with z(T — T c ) = 4.7 ± 0.1, where 
the error is estimated by shifting the data until they visibly 
do not agree. The dashed line is a guide to the eye. 
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FIG. 19: (Color online) Spin-glass susceptibility from equi- 
librium measurements x(L) at T = 0.434 in three dimen- 
sions. x(L*) is the susceptibility determined from the off- 
equilibrium simulations with L* = At 1/Iz , A = 1.27, and 
z « 4.67. We choose L* so that the data for x{L) and x(^*) 
fall on a straight line allowing us to determine the dynamical 
critical effective exponent z(T) as a function of temperature. 
When L*(t) approaches the sample size (here L = 16), x(^*) 
necessarily saturates. The dashed line is a guide to the eye. 
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FIG. 21: (Color online) Data for the autocorrelation function 
q(t) in three dimensions as a function of time (measured in 
Monte Carlo steps). The dashed line is a guide to the eye to 
emphasize the asymptotic slope of q(t) in a log-log plot. 



timated rj*(T) are very similar over a range of temper- 
atures, as shown in Fig. ^| Hence it is not possible to 
obtain an independent measurement of T c using the con- 
sistency criterion outlined above (see Fig.l21|). 
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FIG. 22: (Color online) Scaling plot of the data for &l/2£ l 
vs £l/L in two dimensions. The deviations at large £,l/L 
suggest corrections to scaling. If T c is finite, the curves must 
go through £ 2 l/2£l = 1 at T = T c . The data only tend 
asymptotically to 1 proving that T > T c . 
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FIG. 23: (Color online) Examples of excellent dynamic scal- 
ing well into the paramagnetic region in two dimensions. 
Spin-glass susceptibility from equilibrium measurements x(£) 
(triangles) and x(L*) from scaled dynamical measurements 
(squares) at two temperatures: T — 0.130 (upper set of 
curves) and T = 0.309 (lower set of curves). 



VII. TWO DIMENSIONS 



The gauge glass in two dimensions has been studied 
extensively in Refs.0 and I 111 In particular it was found 
that for this model T c = 0, inferred from the values of 
critical exponents and several tests made. In this section 
we present results for off-equilibrium simulations, as well 
as another test which shows that the system orders at 
zero temperature in two dimensions. 



A. Correlation length 

The GG in dimension two has been extensively 
studiedA2^i^2iM44i^ Measurements on various observ- 
ablcs, in particular the correlation length, have estab- 
lished that the ordering temperature is either zero or 
much lower than our lowest measuring temperature, T = 
0.13iii In Fig. we show correlation length data pre- 
sented in the form of plots of the ratio £(2Z>)/2£(L) = 
^2l/2^l against £(L)/L. In the absence of corrections to 
scaling all points should fall on a single scaling curve, 
which for a finite ordering temperature should pass 
through £(2L)/2f(L) = 1 at T = T c [here £(L)/L is 
independent of L] ^ For ordering which takes place only 
at zero temperature the ratio should tend asymptotically 
to a value which can be less than 1. The data show this 
form of behavior, demonstrating conclusively that for this 
system T c = 0. 



B. Nonequilibrium scaling 



Measurements of the equilibrium susceptibility x(-^) 
for L = 4 to 16 (or to 24 above T = 0.20) and of 
the nonequilibrium susceptibility x(t) f° r samples of size 
L — 64 are shown in Fig. [23 As for the other dimensions, 
the nonequilibrium data are scaled using L* — At x l z 
with A and z chosen at each temperature so that the 
two sets of susceptibility data scale together. The scal- 
ing is particularly satisfactory because a wide range of 
L could be used for the equilibrium measurements, and 
because the nonequilibrium sample sizes have been cho- 
sen such that for all measurements L* is much smaller 
than the sample size, so there are no saturation effects 
at long times. This demonstrates once again, but more 
clearly than for higher dimensions, that a dynamical ex- 
ponent z(T) can be defined relating the annealing time t 
to an effective length L*(t) for temperatures well above 
any critical temperature. The scaling holds to within 
the present statistical accuracy for the whole range of 
L from L — 4 to L — 24, or alternatively for anneal 
times from 2 MCS to over 16000 MCS. It appears that 
the dynamic scaling concept is ubiquitous and is not just 
valid at critical points. Finally, Fig. [2] shows data for 
the autocorrelation function q(t). One can clearly see 
that the data, presented in a log-log plot, are strongly 
curved suggesting that the standard scaling functions for 
a finite-temperature transition do not apply. 
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FIG. 24: (Color online) Data for the autocorrelation func- 
tion q(t) in two dimensions as a function of time (measured 
in Monte Carlo steps). Here L = 64 and T = 0.13. The 
curvature of the log-log plot demonstrates yet again that any 
freezing is well below this temperature. 



VIII. COMPARISON OF TECHNIQUES 

We can draw comparisons between the different tech- 
niques for estimating T c and the critical exponents. It 
should be noted that equilibrium data for the different 
observables studied are taken within the same runs which 
allows a comparison of the relative statistical precision of 
the different measurements for exactly the same compu- 
tational effort. As an example, for measurements made 
at L — 8 at T = 0.461 in dimension three and with 
the simulation parameters listed in Table [HJ the spin- 
glass susceptibility x 1S accurate to 0.5%, the correlation 
length £(Z) to 0.7%, the current I lms to 1.6% (and the 
internal energy to 0.02%) with the quoted errors purely 
statistical. 

The major difficulty in the interpretation of the data 
is not statistics but corrections to scaling. In the gauge 
glass the current method turns out to be relatively insen- 
sitive to corrections to finite-size scaling and is therefore 
already reliable for relatively small samples. Averaging 
over large numbers of samples is, however, essential for 
these measurements because of strong intrinsic sample to 
sample fluctuations. In ISGs measurements of the critical 
behavior of the domain-wall stiffness (see, for instance, 
Ref. l46|) could well provide an equally efficient and reli- 
able method to estimate T c , applicable quite generally. 

Direct analysis of the spin-glass susceptibility x{L,T) 
on its own is also reliable and efficient, particularly as the 
statistical precision of this parameter is high. At T = T c 
there is pure power-law behavior, x(L,T)/L 2 ~ L~ V ^ T \ 
The onset of deviations from the pure power-law behavior 



below and above the ordering temperature provides a 
clear signature for T c . 

In dimension four, even with the restricted range of L 
available, there is a downward deviation from the pure 
power law as soon as T > T c , and an upward deviation 
from the pure power law for temperatures T < T c . Cor- 
rections appear to be negligible as the straight line fit at 
T c is excellent even for L — 2. 

In dimension three the onset of deviation from a pure 
power law to a downward curvature at temperatures 
above T ~ 0.46 gives a clear indication of T c but the 
low-temperature upturn is weak. Corrections are present 
as there is a weak curvature at small L for all tempera- 
tures. If the plausible assumption is made that the lead- 
ing correction to scaling is the leading analytic term, a 
precise estimate for T c is obtained from fits using all the 
x(L,T) data over the whole range of L. This estimate is 
consistent with the value from the current measurements. 



The susceptibility method (see Sees. IV Bl and rVTBf) 

becomes progressively more precise when data for larger 
samples are available. In the present case, data up to L — 
5 only are sufficient for dimension four, but in dimension 
three it would have been very helpful to have had data up 
to L = 12; with such data in hand it would be possible 
to pinpoint the transition temperature to even higher 
precision from susceptibility measurements alone. 

The technique involving the consistency of observables 
deduced from nonequilibrium scaling and equilibrium 
susceptibility gives an estimate for T c which has been 
found to be in excellent agreement with those obtained 
from the other methods for dimension four. This de- 
spite the presence of short-time corrections to scaling 
in q(t) that were unexpected. In dimension three con- 
sistency between alternative methods of estimating pa- 
rameters holds over a range of temperatures rather than 
at a unique temperature which could be identified with 
T c . This accidental consistency has not been observed in 
other systems and remains to be understood, but means 
that the technique in this particular case cannot be used 
to estimate T c with any precision. The nonequilibrium 
measurements are relatively economical in computational 
resources as full thermodynamic equilibration is not re- 
quired. 

Estimates of T c from the intersections of correlation 
length ratios £(L)/L turn out to be misleading in dimen- 
sion four because of strong intrinsic corrections to finite- 
size scaling for the small or moderate system sizes to 
which computational resources generally limit measure- 
ments in spin-glass systems at high dimensions. This pa- 
rameter thus gives a qualitative indication that spin-glass 
ordering is occurring but it is not reliable for extracting 
precise values of the ordering temperature from small L 
data even in a situation where other measurements can 
give consistent and satisfactory estimates for an equiva- 
lent range of L. It can be noted that while corrections to 
finite-size scaling in measurements of x(L) appear to be 
negligible for the GG in dimension four, the corrections 
for are still very strong at size L — 5. Even in 
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TABLE V: Critical temperature T c and critical exponents 
for the GG in different space dimensions d and from different 
references. In the table KC represents results from this work, 
whereas RY represents Ref . by Reger and Young, OY Ref . [fj 
by Olson and Young, and K Ref. [ll] by Katzgraber. 



d 


Reference 


T c 


V 


V 


Z LO 


4 


RY 


0.96(1) 




0.70(15) 




4 


KC 


0.89(1) 


-0.74(3) 


0.70(1) 


4.50(5) 


3 


OY 


0.47(3) 


-0.47(7) 


1.39(20) 


4.2(6) 


3 


KC 


0.460(15) 


-0.47(2) 


1.39(5) 


4.7(1) 2.5(5) 


2 


K 








2.56(20) 


oo 



TABLE VI: Critical temperature T c and critical exponents 
for the ISG in different space dimensions d with Gaussian 
interactions. In the table PRR represents R ef. [491 by Parisi, 
Ricci-Tersenghi and Ruiz-Lorenzo, MPR Ref.|5Cjby Marinari, 
Parisi, and Ruiz-Lorenzo, CEA Ref. l5ll by Campbell et at, 
and BY Ref. by Bhatt and Young. 



d 


Reference 


T c 


V 


V 


z 


4 


PRR 


1.80(1) 


-0.35(5) 


1.0(1) 




4 


CEA 


1.78(1) 


-0.44(2) 


1.08(10) 


4.9(4) 


3 


MPR 


0.95(4) 


-0.36(6) 


2.00(15) 




3 


CEA 


0.92(2) 


-0.42(3) 


1.65(5) 


6.45(10) 


2 


BY 








3.63(10) 





dimension three where the corrections are weaker, this 
observable remains a poor tool for estimating the tran- 
sition temperature. Strong correction effects in £(L)/L 
can also be seen in data reported for XY spin glasses^ 
and have been found in certain ISGs. Therefore the cor- 
rection effect appears to be generic and this technique 
should be applied only with caution. 

Finally, estimates from the Binder parameter are in- 
appropriate for the GG as the curves for different system 
sizes do not intersect for the range of L over which mea- 
surements are carried out. It is possible that there are 
strong finite-size corrections for this parameter so that 
intersections would only be seen for much larger L. 

Once T c has been estimated, the susceptibility mea- 
surements at that temperature give immediately the cor- 
responding value of the critical exponent r\. The preci- 
sion for this parameter is limited entirely by the accuracy 
with which T c is known. The exponent v is estimated 
from scaling plots of current data and from scaling plots 
of susceptibility data. Once again precision is limited by 
the accuracy with which T c has been determined. The 
dynamical exponent z can be measured accurately from 
a comparison of equilibrium and nonequilibrium suscep- 
tibilities. The data show that z(T) can be estimated 
operationally over a wide range of T including temper- 
atures well above the ordering temperature [z{T) varies 
continuously through the ordering temperature] . 

In the light of this extensive analysis, it would seem 
appropriate to critically re-assess the estimates for the 
critical exponents in the ISG family. 

IX. CRITICAL PARAMETERS AS A 
FUNCTION OF SPACE DIMENSION 

The critical parameters for the GG obtained from this 
and earlier work are listed in Table El 

In both dimensions three and four the GG ordering 
temperatures T c are roughly half of those for the Gaus- 
sian ISG. By interpolation between the zero-temperature 
current stiffness exponents 9(d) in dimensions 1, 2, and 
we can estimate the lower critical dimension for the 
GG as the dimension at which 9(d) passes through zero: 
died ~ 2.5, very similar to that of the ISG systems with 



continuous interaction distributions^^ (see Table IVII 

and Refs.mHUll). 

The corrections to scaling in the GG follow just the 
same pattern as in the bimodal (± J) ISGs. In the Gaus- 
sian ISG the corrections to scaling for x(L) at T c appear 
very weak and so unmcasurablc. In dimension four for 
both cases the correction to scaling for the susceptibil- 
ity x(L) is so weak as to be unobservable even down 
to L = 2. In dimension three there is a clear correc- 
tion to scaling for x(L) which can be fitted satisfac- 
torily to the form ~ (1 + AL~" rfT ). w off in the ISG 
case has been interpreted as the leading irrelevant op- 
erator correction^ and this could also be the case for 
the three-dimensional GG. However, it seems more likely 
that the correction is dominated by a "lattice artifact" 
giving cj c ff = 2 — 77, which has no relation to the leading 
irrelevant operatorii2ii£ For both the three-dimensional 
GG and the three-dimensional bimodal ISG, the value of 
the apparent correction to scaling exponent is compati- 
ble with this interpretation, which would imply that the 
leading irrelevant operator term has an exponent value 
uj > 2 - 77. 

Finally, we can give an overview of the critical expo- 
nents in the GG and ISG families as a function of space 
dimension. It should be noted that there are some con- 
straints that apply to both families: at the upper critical 
dimension d = 6, v = 1/2, r\ = 0, and z = 4. If T c = 0, 
v = —1/9, where 9 is the stiffness exponent. There- 
fore, at the lower critical dimension where the stiffness 
exponent 9 = 0, v is infinite. For systems with a nondc- 
generate ground state (such as the GG or the ISG with 
Gaussian interactions, but not the ISG with bimodal in- 
teractions), when T = there is no decay of the correla- 
tion function G(r) with r. This necessarily implies that 
when T c = 0, rj(d) = 2 — d. [For systems with degenerate 
ground states rj(d) > 2 — d]. 

The leading terms in the RG e-expansion for the GG 
are u(d) = 1/2 + be/24 and rj(d) = -e/6. For the ISG 
u(d) = 1/2 + 5e/12 and rj(d) = It can be seen 

in Figs. 1251 and !26l respectively, that in the GG as well as 
for the ISG the leading e-expansion term gives the right 
sign and a qualitative indication of the strength of the 
variations of the exponents just below the upper criti- 
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FIG. 25: (Color online) Critical exponent v as a function of 
space dimension for the ISG with Gaussian bonds and the 
GG (data taken from Tables IVT1 and IV! respectively). The 
shaded region denotes the area where one expects the lower 
critical dimension. The solid [dashed] line represents the RG 
e-expansion estimate for the GG [ISG], as shown in the text, 
and the dotted lines are guides to the eye. At the lower critical 
dimension one expects u(d — » di c d) — * oo. The data for the 
ISG and GG support this behavior. 



cal dimension. The leading terms in the expansion for 
■>] cannot, however, be used to predict that the rj{d) val- 
ues are much more negative in the GG than in the ISG. 
The observed exponents deviate strongly from the curves 
calculated to third order for the ISG where the second 
and third order expansion terms are very large. This de- 
viation is smaller for the GG. The leading term in the 
e-expansion for the exponent v works better, with v(d) 
being higher for the GG than for the ISG, in agreement 
with the relative strengths of the leading terms (includ- 
ing higher order terms would entirely destroy the agree- 
ment). This is all in striking contrast to the canonical 
ferromagnets (Ising, XY, or Heisenberg) without disor- 
der, where the third order e-expansion correction gives 
excellent predictions for 77(d) an v(d). 

If we look at the overall behavior of the exponents over 
the whole range of d, the shape of the v(d) curves is fairly 
similar for the GG and Gaussian ISG families, although 
the divergence at the lower critical dimension di cc i « 2.5 
is distinctly narrower in the GG case. The shapes of the 
77(d) curves are on the other hand dramatically different. 
The leading terms in the e-expansions agree poorly with 
the data even for d = 4, with an inversion of the mea- 
sured positions of 77(d) as compared with the expansion 
predictions. For the Gaussian ISG family 77(d) decreases 
as the dimension drops from the upper critical dimen- 
sion to the lower critical dimension. For the GG there 
is a deep minimum in 77(d) somewhere in the region of 




d 



FIG. 26: (Color online) Critical exponent rj as a function 
of space dimension. The diagonal dotted line represents the 
physical limit line 77(d) = 2 — d, whereas the horizontal dotted 
line represents 77 = 0. The solid and dashed lines represent the 
e-expansion estimates for 77 for the GG and ISG, respectively. 
For all d < di c d ~ 2.5 77(d) must join the limit line. The 
shaded region denotes the area where one expects the lower 
critical dimension. 

d = 4 (data at d = 5 would be needed to pin down 
the position of the minimum). In spin glasses, simply 
going from Ising to vector spins changes this exponent 
considerably, whereas in three-dimensional ferromagnets 
r] is very small and practically independent of the type 
of spin symmetry. 

Finally, in Fig.[^|the critical dynamical exponent z(d) 
can be compared near d — 6 with the van Hove approx- 
imation z = 2(2 - or for the ISG case with the 
first order e-expansion z = 2(2 — rj)/(l + e/4)i£& The van 
Hove expression can provide a qualitative indication for 
z(d) once the numerical values of 77(d) are known, but if 
the directly measured 77(d) values are used rather than 
the epsilon expansion estimates, this expression would 
predict higher z values at each d for the GG compared 
with the ISG values, which is not what is observed. z(d) 
must diverge at the lower critical dimension, while curi- 
ously the GG z{d) values only increase slightly between 
d = 4 and d = 3. For the Gaussian ISG, z(d) shows an 
indication of divergence as the space dimension drops. 



X. CONCLUDING REMARKS 

We have studied the critical parameters of the GG, 
in dimensions two, three and four. The data confirm 
that T c is zero in two dimensions. It appears that the 
most reliable and accurate methods for estimating critical 
temperature and exponent values are direct susceptibil- 
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FIG. 27: (Color online) Critical exponent z as a function of 
space dimension. The horizontal dotted line marks the mean- 
field asymptotic value z — 4, whereas the shaded region de- 
notes the area where one expects the lower critical dimension. 
The solid (dashed) curve shows the asymptotic behavior as 
predicted from the van Hove approximation (Ref. 55) for the 
GG (ISG). One expects that z(d -> di cd ) -> oo. The ISG data 
show this behavior whereas the GG data are inconclusive. 

ity measurements, effective stiffness measurements (cur- 
rent measurements in the GG case), and a combination 
of nonequilibrium and susceptibility measurements. In 
each case corrections to scaling are either so small as to 
be negligible or can be taken into account. Binder pa- 
rameter measurements are inappropriate as the Binder 
parameter curves do not intersect, at least for the sys- 
tem sizes used, and correlation length measurements are 
not always reliable because there can be strong deviations 
from finite-size scaling. 

From these results a general strategy for obtaining 
(high-precision) critical parameters numerically in any 
glassy system can be sketched out. It is essential to 
be able to rely on high-quality equilibrium susceptibility 
data over a wide range of sizes and to as large a size as 
computational resources allow. Careful analysis of these 
data allowing for the leading correction to scaling term 
should provide reliable estimates of the ordering temper- 
ature and of the equilibrium critical exponents. Further 
measurements of the effective stiffness as a function of 
temperature together with nonequilibrium behavior can 



then confirm the value of the ordering temperature and 
give further information on other exponents including 
the critical dynamical exponent z. Parameters involv- 
ing ratios (such as the Binder parameter or the correla- 
tion length) seem frequently to be biased, at small and 
moderate sizes, by complicated corrections to scaling and 
should be treated with caution. 

If we compare with other systems having d = 6 as 
upper critical dimension, ISGs and percolation, the ex- 
ponent v(d) always evolves regularly, diverging when the 
dimension reaches the lower critical dimension. The ex- 
ponent rj{d) on the other hand changes dramatically in 
form from one family to the next. For the Gaussian 
ISG rj{d) grows smoothly more negative as the dimen- 
sion drops from the upper critical dimension to the lower 
critical dimension. For the percolation systems 77(d) has a 
weak minimum around d — 3 before becoming positive at 
d = 2. 57 The present results show that for the GG there 
is a deep minimum in 77(d) near d = 4. This is in stark 
contrast to the situation in the canonical Ising, XY, or 
Heisenberg systems with no disorder where rj(d) hardly 
changes at all when the degrees of freedom of the spin are 
modified. In addition, for the canonical systems the RG 
e-expansion to third order gives excellent predictions for 
the exponents at dimensions well below the upper crit- 
ical dimension, while for the ISG and GG systems the 
e-expansion to the same order^^4 gives poor predictions 
for dimensions quite close to the upper critical dimension 
d = 6. 

We have thus proposed a road plan for determining 
critical exponents reliably and accurately from simula- 
tions; however, the results confirm once again that the 
canonical RG theory lacks essential ingredients when ap- 
plied to glassy systems. The physical significance of 
the exponent values obtained, however good they are, 
must remain obscure until the appropriate theoretical 
approach going beyond the traditional RG is found for 
interpreting them. 
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